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摘 要 :为 研究 损伤 阅 值 对 存在 天 然 缺 陷 软 岩 损伤 力学 特性 的 影响 ,基于 损伤 阅 值 理论 ,结合 修正 
的 Lade-Duncan 强度 准则 与 弹性 胡 克 定律 ,推导 出 软 岩 损伤 国 值 表 达 式 ,利用 Weibull 分 布 理论 及 
Lemaitre 应 变 等 价 理论 构建 软 岩 统计 损伤 本 构 模 型 ,同时 结合 三 轴 压 缩 试 验 对 其 开展 验证 。 研 究 结 
果 表 明 : 基 于 修正 的 Lade-Duncan 强度 准则 与 弹性 胡 克 定律 求解 出 的 软 岩 损伤 阅 值 过 于 保守 ,同时 
天 然 缺 陷 对 岩石 的 轴 向 应 变 影 响 较 大 ;以 岩石 中 天 然 缺 陷 为 修正 变量 ,对 损伤 阔 值 进行 修正 后 , 软 
兰 的 损伤 变化 规律 破坏 过 程 描述 及 应 力 - 应 变 演 化 更 贴 合 实际 ,验证 了 本 研究 软 岩 损伤 本 构 模 型 


的 科学 性 、 合 理性 。 
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Abstract: lo study the influence of damage threshold on the damage mechanical properties of soft rock 
with natural defects, based on the damage threshold theory, combined with the modified Lade-Duncan 
strength criterion and the elastic Hooke's law, we derived the expression of the damage threshold of soft 
rock. The statistical damage constitutive model of soft rock is constructed by using Weibull distribution the- 
ory and Lemaitre strain equivalence theory, and verified by a triaxial compression test. The results show 
that the damage threshold of soft rock based on the modified Lade-Duncan strength criterion and elastic 
Hooke's law is too conservative ,and the natural defects have a great influence on the axial strain of rock. 
Taking the natural defects in rock as modified variables ,the damage variation law , failure process descrip- 


tion , and stress-strain evolution of soft rock are more practical after modifying the damage threshold , which 
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verifies the scientificity and rationality of the damage constitutive model of soft rock in this paper. 


Key words:Lade-Duncan strength criterion ; damage threshold ; statistical damage model;rock mechanics; 


three parameter Weibull distribution 


岩石 是 一 种 非 均 质 的 天 然 地 质 材料 ,页 岩 、 砂 岩 
等 含油 气 岩 石 是 典型 的 代表 ,其 内 部 存在 很 多 分 布 
不 规律 的 微 裂 颖 和 空隙 ,由 于 微观 缺陷 使 得 岩石 的 
损伤 破坏 机 制 变 得 很 复杂 ,在 岩石 发 生 损 伤 的 过 程 
中 ,其 内 部 微 有 裂 颖 的 宽度 与 其 长 度 存在 复杂 的 关 
RU ,孔隙 压力 的 增加 也 会 造成 有 效应 力 系数 的 减 
小 。 岩 石 微观 缺陷 可 视 作 为 损伤 的 实质 性 表现 ， 
针对 于 岩石 本 构 模 型 的 研究 ,是 岩石 力学 研究 领域 


规律 及 破坏 特征 进行 研究 ,并 建立 岩石 损伤 模型 ; 梁 
明 纯 等 "考虑 岩石 在 发 生 损伤 破坏 过 程 中 所 表现 
出 来 的 剪 胀 扩容 特性 ,基于 损伤 力学 建立 损伤 模型 。 

学 者 从 不 同 的 角度 建立 的 岩石 损伤 本 构 模 型 ， 
在 一 定 程度 上 推动 了 岩石 力学 的 发 展 。 但 是 随 着 研 
究 的 不 断 深 入 ,学 者 对 岩石 的 统计 损伤 模型 进行 了 
修正 及 优化 ,使 其 能 够 更 好 地 描述 岩石 损伤 演化 规 
律 以 及 岩石 的 力学 特性 。 针 对 岩石 损伤 模型 的 修正 


中 的 一 大 重点 ,而 现 有 针对 于 岩石 损伤 破坏 本 构 模 
型 的 研究 ,多数 是 基于 统计 损伤 理论 。 

关于 岩石 的 损伤 模型 ,前 人 从 不 同 的 角度 建立 
了 岩石 的 损伤 本 构 模 型 。 关 顺 等 所 分析 了 现 有 的 分 
数 阶 流 变 本 构 模型 在 描述 岩石 损伤 破坏 时 的 不 足 ， 
在 分 数 阶 微 积分 理论 和 统计 损伤 力学 理论 基础 上 ， 
建立 了 简洁 且 准 确 的 描述 泥岩 蠕 变 的 新 型 损伤 本 构 
模型 ,并 对 其 进行 了 验证 ; 张 超 等 趾 从 岩石 微 元 破坏 
的 角度 分 对 岩石 破坏 进行 了 研究 ,并 基于 微观 损伤 
理论 提出 了 已 经 破坏 的 微 元 体会 对 岩石 引起 非 均 匀 
性 的 函数 假设 ,在 该 假设 的 基础 上 提出 了 岩石 的 损 
伤 模 型 ; 刘 文 博 5 gue" ugue" Wee 
坏 过 程 中 的 能 量 消耗 角度 出 发 ,探讨 岩石 损伤 演化 
规律 及 破坏 变形 特征 ,并 提出 了 能 够 较 好 地 描述 岩 
石 损伤 破坏 过 程 及 其 力学 性 能 的 损伤 本 构 模 型 ; 
SUN 等 中 基于 压 裂 曲线 的 原 位 脆性 评价 方法 ,推导 
了 页 岩 在 现场 应 力 和 温度 下 的 原 位 脆性 ,并 验证 了 
可 以 通过 观察 压 裂 曲线 的 波动 来 了 解 岩 石 的 压 裂 特 
征 等 损伤 状况 ;MA 等 中 利用 X 射线 CT( 计 算 机 断 
层 扫描 ) 扫 描 技 术 和 损伤 力学 ,提出 了 一 种 识别 水 合 
页 岩 损伤 特征 的 新 型 表征 方法 ; YU 等 中 基于 损伤 
力学 建立 了 渗透 压 下 岩石 损伤 本 构 模 型 ;LIU 807 
定义 了 一 种 新 颖 的 双重 损伤 热力 学 模型 ,并 提出 了 
双重 损伤 本 构 模 型 ;李兵 大 '"" 、 欧 雪 峰 等 后 从 冲击 
动力 学 的 角度 对 岩石 的 损伤 演化 规律 进行 了 研究 ， 
并 建立 了 能 够 描述 岩石 应 力 -应 变 的 损伤 本 构 模 型 ; 
王 东 等 “5 基于 岩石 的 破坏 类 型 建立 其 损伤 模型 ， 
并 通过 岩石 的 单 轴 及 三 轴 压 缩 试验 进行 验证 。 刘 德 
EUS 余 洁 等 .1 对 于 不 同 围 压 下 岩石 的 损伤 演化 


主要 从 以 下 几 个 方面 进行 :中 针对 岩石 的 残余 强度 
进行 修正 。XU! zy axo 通过 原 
有 的 损伤 模型 进行 分 析 , 发 现 其 在 峰 后 以 及 岩石 残 
余 强 度 的 描述 不 足 ,引入 相关 的 修正 系数 对 岩石 峰 
后 阶段 及 其 完全 破坏 后 的 残余 强度 进行 了 修正 ; 
@ 从 岩石 微 元 体 的 损伤 破坏 对 岩石 整体 破坏 的 影响 
进行 修正 , 张 超 等 所 考虑 岩石 微 元 体 发 生 损伤 时 ,对 
于 岩石 整体 的 影响 提出 一 个 修正 系数 ;从 假设 岩 
石 微 元 破坏 所 服从 的 强度 准则 出 发 进行 修正 ， 
TIAN 7 、 田 振 元 等 ”通过 引入 修正 后 的 强度 准则 
对 岩石 损伤 模型 进行 修正 ;由 针对 岩石 的 不 同 赋 存 
条 件 进行 修正 , 视 凯 等 .站 针对 岩石 在 不 同 的 赋 存 条 
件 下 的 岩石 本 构 模 型 的 研究 ,提出 一 个 参数 对 损伤 
模型 进行 修正 ;@LI 等 ”加 入 反映 应 力 下 降 速 率 的 
修正 因子 ,建立 了 考虑 应 力 下 降 特性 的 深 脆 性 岩石 
损伤 本 构 模 型 ;@QU 等 “引入 损伤 校正 系数 ,对 化 
学 冻 融 砂岩 在 荷载 作用 下 的 总 损伤 进行 修正 。 

关于 岩石 的 损伤 本 构 模 型 ,以 及 对 其 进行 相关 
修正 的 研究 ,国内 外 的 研究 学 者 取得 了 丰富 的 成 果 。 
基于 修正 Lade-Duncan 强度 准则 来 确定 软 岩 损伤 阔 
值 ,并 建立 损伤 本 构 模 型 的 研究 成 果 较 少 。 本 研究 
将 采用 基于 Lemaitre 应 变 等 价 假设 理论 ,考虑 损伤 
国 值 对 软 贿 力学 特性 的 影响 ,结合 修正 Lade-Duncan 
强度 准则 和 弹性 胡 克 定律 推导 损伤 阔 值 的 求解 表达 
式 , 同 时 结合 三 轴 压 缩 试验 对 求解 表达 式 进行 验证 ， 
运用 数理 统计 分 析 中 的 多 元 线性 回归 拟 合 的 方法 求 
解 出 修正 系数 表达 式 , 认 为 软 岩 微 元 的 损伤 破坏 服 
从 三 参数 的 Weibull 分 布 ,以 此 建立 软 岩 的 损伤 本 构 
模型 ,引用 三 轴 压 缩 试验 对 其 进行 验证 。 
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1 试 统计 损伤 本 构 模 型 


由 于 岩石 的 非 均 匀 性 ,为 了 能 够 更 好 地 描述 岩 
石 的 破坏 过 程 及 其 力学 特性 ,所 以 当前 针对 岩石 的 
损伤 本 构 模 型 多 数 是 基于 概率 统计 理论 和 损伤 力学 
理论 而 开展 研究 的 。 


1.1 损伤 演化 函数 


岩石 发 生 损伤 的 过 程 可 看 作为 岩石 微 元 逐渐 发 
生 损 伤 并 积累 的 过 程 , 当 作 用 荷载 达到 某 一 水 平时 ， 
将 岩石 整体 划分 为 NV 个 微 元 , 设 已 经 发 生 破坏 的 微 
元 数 为 N,, 定 义 统计 损伤 变量 D 为 发 生 破坏 的 微 元 
数 与 总 微 元 数 比 值 的 表达 式 , 即 
pan a) 
N, 
在 岩石 加 载 的 弹性 阶段 ,认为 其 还 未 发 生 损伤 ， 
此 时 损伤 变量 D =0。 随 着 荷载 的 逐渐 增加 ,应 力也 
随 着 增 大 ,此 时 其 内 部 的 微 裂缝 开始 逐渐 扩展 ,产生 
新 的 裂 颖 ,岩石 也 开始 发 生 损伤 , 即 岩 石 微 元 发 生 损 
伤 。 认 为 岩石 微 元 强度 F IRA Weibull 概率 分 布 ， 
任意 应 变 水 平 区 间 [ c, sy + de ] 中 发 生 破坏 的 微 
元 数目 为 


dN, = Np(ei)de (2) 

由 于 岩石 损伤 破坏 机 制 的 复杂 性 ,为 了 考虑 损 

伤 阔 值 对 岩石 力学 特性 的 影响 ,采用 三 参数 的 

Weibull 概率 分 布 , 其 比 常规 的 Weibull 分 布 多 了 一 

个 损伤 阔 值 参数 sis 。 当 岩石 的 轴 向 应 变 se, KEH 

HERE ca 时 , 其 损伤 量 D 的 概率 分 布 函 数 表 达 

式 为 
m fer E £j — u 

p(ei) = A exp| ? [5 | 


20 
(3) 
式 中 ,m 和 es 分 别 为 概率 分 布 函数 的 形状 参数 和 尺 
设 岩 石 随 围 压 加 载 到 某 一 应 力 水 平时 ,其 已 经 
发 生 损 伤 的 微 元 数目 NS 


Ni = | Np( ei) de, (4) 
则 损伤 变量 为 
0,el S en 
D= E m (5) 
1 — exp| 一 [f 09) IET > £u 
£o 
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1.2 Lade-Duncan 强度 准则 


现 有 常用 的 岩石 强度 准则 ,Morh-Coulomb 准则 
表达 式 人 简单 方便 ,但 其 未 考虑 中 主 应 力 的 影响 ;而 
Drucker-Prager 准则 是 对 Morh-Coulomb 准则 的 近似 ， 
考虑 了 中 主 应 力 的 影响 ,但 其 确定 的 岩石 微 元 强度 
的 合理 性 有 限 ; Hoek-Brown 准则 是 根据 现场 经 验 而 
建立 的 ,其 随意 性 较 大 ; 而 Lade-Duncan 准则 是 针对 
无 黏 性 土建 立 的 ,经 过 修正 可 以 运用 于 c- o 材料 , 合 
理 地 考虑 了 中 主 应 力 的 影响 ,可 以 较 好 地 描述 岩石 
的 力学 行为 。 岩 石 微 元 破坏 强度 准则 服从 Lade- 
Duncan 强度 准则 , 即 
n (3 - sing)? 

L TO C] + TE sing)? 


(6) 


pK 为 岩石 的 材料 常数 ; e C) NAERA; T. 
分 别 为 应 力 第 一 、 三 次 不 变量 ,表达 式 为 


p = Oil +0, +0, 


(7) 


l, = 01003 

式 中 , cl o o (MPa) 分 别 为 大 、 中 、 小 主 应 力 。 

该 强度 准则 是 由 Lade 和 Duncan 对 大 量 的 沙土 
三 轴 压 缩 试验 结果 的 分 析 总 结 得 出 的 关于 沙土 破坏 
强度 的 理论 表达 式 , 其 主要 适用 于 无 黏 性 土 。 为 了 
使 其 能 够 描述 c- o 材料 的 力学 特性 ,扩展 其 适用 范 
围 ,EWY 等 ”通过 引入 黏 结 应 力 对 其 进行 了 修正 。 

gy = ccotp (8) 

式 中 ,ce HAARR o 

修正 后 的 Lade-Duncan 强度 准则 表达 式 为 


cra. (3 - sing)? (9) 
L CI + sing) C1 - sing)? 
y = Ol +0, +0; +309 (10) 
L -6g, +0) T * a9) (o4 * o) 


WE X C9) 和 式 (10) ,可 得 出 岩石 材料 的 Lade- 
Duncan 强度 准则 表达 式 为 
Oil tO, 1 4 3 3 2 
人 
(11) 
ol = 90, + (nm - Dos (12) 


1.3. 损伤 阔 值 确定 

宕 石 在 加 载 过 程 中 , 当 洗 石 的 轴 疝 应 变 未 达到 
损伤 浆 值 时 ,此 时 岩石 微 元 处 于 弹性 阶段 且 满 足 弹 
性 胡 克 定律 。 
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el = [o, - vg, 4 g,) VE (13) ”应 的 轴 向 应 变 , 这 与 实际 试验 是 不 符合 的 ; 当 p = 
在 常规 的 三 轴 试 验 中 go,= 0, ,得 26.5" 时 ,此 时 求 出 来 的 损伤 阔 值 也 是 过 于 保守 的 ， 

& = [ol - 2vug,V/E (14) ”但 是 未 出 现 异常 数据 。 


利用 岩石 在 损伤 国 值 处 的 特殊 性 , 当 岩 石 的 轴 
向 应 变 处 于 损伤 阔 值 时 ,岩石 处 于 即将 开始 发 生 损 
伤 的 临界 状态 ,此 时 的 岩石 仍 处 于 弹性 阶段 , 联 立 弹 
性 胡 克 定律 和 修正 的 Lade-Duncan 强度 准则 ,推导 
出 损伤 阔 值 的 求解 表达 式 , 即 
em = [go + (q- Do, -2vao,VVE (15) 
将 式 (15) 代 入 式 (5), 即 可 得 到 损伤 变量 D 与 
轴 疝 应 变 e, 的 关系 表达 式 。 引 用 软 岩 的 试验 数据 
对 式 (15 ) 进行 验证 ,材料 参数 和 力学 参数 分 别 如 
表 1 和 表 2 所 示 。 软 宕 1、 软 宕 2 I E UBER LK 
特点 是 质量 轻 、 孔 洞 多 、 硬 度 不 高 。 其 力学 行为 具有 
弹 塑 性 应变 硬化 .应变 软化 和 时 间 依 赖 性 。 无 侧 限 
抗 压 强度 为 1 ~20 MPa, 其 力学 行为 介 于 土壤 和 岩石 
之 间 。 如 图 1(a) 所 示 , 凝 灰 岩 在 相对 低 倍 下 呈 碎 悄 
集合 状 和 团 簇 状 ;如 图 1(b) 所 示 , 薄 层 被 放大 后 , 主 
体 呈 弯曲 片 层 状 结构 。 
表 1 软 岩 的 材料 参数 


Tab.1 Material parameters of soft rock 


岩石 Qt B. p Br a b 


fu 108) — 9.135 0.241 2.657 0.792 0.092 1.085 
fu2U9) — 7.933 0.184 4.211 0.597 1.085 0.092 
R2 软 岩 的 力学 参数 


Tab.2 Mechanical parameters of soft rock 


岩石 E/MPa v C/MPa e/((?) 
软 岩 108 800 0.5 0.56 38.5 
e 209 950 0.5 0.76 26.5 


通过 式 (15 ) 计 算得 到 的 损伤 阐 值 列 于 表 3。 当 


o =38.5° 时 ,运用 修正 Lade-Duncan 强度 准则 求 出 
来 的 损伤 阔 值 过 于 保守 ,接近 于 峰值 应 力 强 度 所 对 
应 的 轴 向 应 变 值 ,并 且 当 围 压 增 大 到 一 定 程 度 时 损 
伤 阔 值 出 现 了 异常 ,损伤 阔 值 大 于 了 峰值 强度 所 对 


Ww 
UA Ls * 


(a) 凝 灰 岩 的 微观 形 貌 (b) 凝 灰 岩 的 片 状 聚 集体 


"i 凝 灰 岩 扫描 电镜 图 


Fig.1 Scanning electron microscopic image of tuff 

针对 以 上 问题 ,通过 对 Lade-Duncan 强度 准则 
以 及 岩石 损伤 破坏 本 构 曲 线 的 特征 分 析 讨 论 ,得 出 
以 下 结论 。 

1) 由 于 Lade-Duncan 强度 准则 一 开始 是 针对 于 
无 黏 性 土 而 提出 的 ,尽管 修正 的 Lade-Duncan 强度 
准则 能 够 用 于 描述 e - o 材料 的 强度 特性 ,但 是 将 该 
强度 准则 运用 到 c - o 材料 时 ,很 大 概率 会 过 高 地 预 
测 其 强度 "” ,导致 运用 该 强度 准则 求 出 来 的 损伤 阔 
值 会 过 于 保守 甚至 出 现 异 常 的 现象 。 根 据 试验 数据 
求 出 的 损伤 阔 值 规律 分 析 而 进行 修正 , 当 内 摩 氛 角 
较 大 时 ,计算 得 到 的 损伤 阔 值 会 距离 峰值 应 力 所 对 
应 的 轴 向 应 变更 近 , 甚 至 在 围 压 较 大 的 时 候 , H 
值 大 于 了 峰值 应 力 强 度 所 对 应 的 轴 向 应 变 ,这 与 实 
际 是 不 相符 的 ,所 以 需要 引入 一 个 关于 内 摩擦 角 o 
的 衰减 系数 4。 

2) 从 岩石 损伤 破坏 本 构 曲 线 特征 分 析 , 岩石 损 
伤 破坏 过 程 可 分 为 空 际 压 密 阶 段 ,弹性 变形 阶段 破 
裂 发 展 阶段 ,应变 软化 阶段 和 塑性 流动 阶段 。 在 围 
压 开 始 施加 的 初始 阶段 ,岩石 正 处 于 空隙 压 密 阶段 ， 
压 密 阶 段 产 生 的 轴 向 应 变 也 是 损伤 阔 值 偏 大 的 一 个 
原因 ,为 削减 这 种 影响 提出 了 一 个 关于 初始 围 压 所 
产生 的 应 变 的 削减 参数 B。 


表 3 修正 前 的 损伤 阅 值 


Tab.3 Correction of the previous damage threshold 


et 软 岩 1081 (gp 238.509) 


软 岩 2.” (p =26. 50°) 


p 


E 


0.5 MPa 1.0 MPa 1.5 MPa 2.0 MPa 0.5 MPa 1.0 MPa 1.5 MPa 2.0 MPa 
损伤 阔 值 ”0.007 665 697 0.010849337 0.014032978 0.017216618 0.003 768875 0.004696984  0.00585712 0.008 177 392 
A 0.01131 0.01177 0.01223 0.012 69 0.010 942 0.01131 0.011 77 0.012 69 

根据 以 上 分 析 讨 论 , 本 研究 针对 以 修正 Lade- fA o MEIJE o, 的 修正 , 即 
Duncan 强度 准则 求 出 来 的 损伤 阔 值 进行 关于 内 摩 e = Ae, - B (16) 


投稿 网 站 :http://cjam. xjtu. edu. en ” 微 信 公众 号 :应 用 力学 学 报 


1388 应 用 力学 学 报 


AP: el 为 修正 后 损伤 阔 值 ;4 = cosp ; B = VE 
E 


中 无 (MPa) 是 岩石 在 空隙 压 密 阶 段 的 一 个 加 权 平 均 
变形 模 量 ,由 试验 数据 可 知 ,其 与 弹性 模 量 存在 一 定 
的 关系 ,下 = CE ,C WE; SIR o, (MPa) 存 在 一 定 的 
非 线性 关系 ,使 用 数理 统计 分 析 中 的 多 元 线性 回归 
拟 合 的 方法 得 如 下 表达 式 。 


第 40 卷 
1) mice 17 拟 合 得 出 
C =1. 077 260, ? - 1.477 640, ° +0. 302 710; + 
0. 302 71 (17) 
2) fica 27 拟 合 得 出 
C =2. 220 270, ° -1.263 960, ° +0. 188 630, + 
0. 048 13 (18) 
由 式 (16) 求 出 的 修正 后 的 损伤 阔 值 见 表 4。 


表 4 修正 后 的 损伤 阐 值 
Tab.4 Modified damage threshold 


Jos 199 ( o =38.5°) 


lot 205 p =26.5°) 


Nem 0.5 MPa 1.0 MPa 1.5 MPa 2.0 MPa 0.5 MPa 1.0 MPa 1.5 MPa 2.0 MPa 
HABE 0.005 74593  0.007912856 0.009 790794 0.010775 567  0.00296431 0.003 595476 0.004359251 0.005 082 74 
Ee 0.01131 0.011 77 0.012 23 0.012 69 0.010 942 0.01131 0.011 77 0.012 69 
由 图 2 可 知 , 围 压 的 存在 可 以 有 效 抑制 岩石 损 KA 2 EOS UE , n] AMA EUR 35105 BO (EZ SEC Wo 
伤 的 产生 。 随 着 围 压 增 大 ,修正 前 的 损伤 阔 值 逐渐 ”实际 的 。 
接近 峰值 应 力 强 度 所 对 应 的 轴 向 应 变 es, 当 围 压 大 nds 
T 1.2 MPa 时 ,修正 前 的 软 岩 1 MAREKET 0012 
2。 ,这 与 实际 试验 数据 不 相符 。 其 一 部 分 原因 是 由 Y 
0010| 一 峰值 应 力 强度 对 应 的 条 向 应 变 e. 


于 Lade-Duncan 强度 准则 会 过 高 预测 岩石 的 强度 ， 
导致 理论 计算 得 到 的 损伤 阔 值 大 于 岩石 实际 的 损伤 
闵 值 ; 男 外 一 部 分 原因 则 是 由 于 压 密 阶 段 产 生 的 轴 
向 应 变 ,导致 理 论 损伤 阔 值 滞后 的 现象 。 修 正 后 的 
损伤 靖 值 相 比 于 修正 前 的 损伤 冰 值 减少 了 25. 0496 
~37.41% ,根据 软 岩 1 模型 验证 可 以 看 出 , 式 (16) 
修正 的 损伤 阔 值 较 符合 试验 数据 。 


0.020 


一 峰值 应 力 强度 对 应 的 轴 疝 应 变 s。 
0.018 
一 修正 前 的 损伤 阔 值 s。 


0.016 
一 一 修正 后 的 损伤 阔 值 si 


0.014 


04 06 08 10 12 Du be 18 20 
围 压 /MPa 


图 2 软 岩 1 损伤 阔 值 修正 对 比 图 


Fig.2 Comparison chart of soft rock 1 damage 


threshold correction 
由 图 3 可 知 , 随 着 围 压 的 增长 , 软 岩 2 Wi I 
值 也 逐渐 增 大 ,其 增长 速率 大 于 峰值 应 力 强 度 所 对 
应 的 轴 向 应 变 e. 的 增长 速率 ,逐渐 接近 于 e. ,但 其 
未 出 现 大 于 e. 的 现象 。 修 正 后 的 损伤 阔 值 相 比 于 
修正 前 的 损伤 立 值 减少 了 21.35% ~37.85% 。 根 据 


一 一 修正 前 的 损伤 阐 值 e,。 


一 一 修正 后 的 损伤 阔 值 sl 


0.5 1.0 l5 2.0 
FE /MPa 


图 3 软 岩 2 mh EEX EET 
Fig.3 Damage threshold correction comparison 


diagram of soft rock 2 
1.4 统计 损伤 本 构 模型 的 建立 


从 微观 的 角度 讨论 分 析 了 光石 的 损伤 机 制 , 岩 
石 的 强度 主要 来 源 于 和 结 、 摩 探 强度 。 岩 石 发 生 损 
伤 的 初始 阶段 ,其 强度 主要 来 源 于 医 结 强度 ; 随 着 傈 
载 的 增加 和 岩石 损伤 的 积累 , 微 裂缝 发 育 并 逐渐 形 
成 局 部 贯通 ,在 力 的 作用 下 将 会 产生 相对 滑动 ,由 于 
贯通 截面 的 不 规则 性 , 必 会 产生 滑动 摩擦 力 , 此 时 洗 
石 的 强度 主要 来 源 于 摩擦 强度 ; 随 着 损伤 的 积累 并 
达到 岩石 的 破坏 强度 值 时 ,岩石 发 生 破坏 ,损伤 变量 
D =1, 处 于 塑性 流动 阶段 。 

由 于 岩石 损伤 机 制 的 复杂 性 ,为 了 简化 模型 参 
数 起 见 , 认 为 残余 剪 切 强度 os 等 于 残余 强度 ,基于 
此 得 到 优化 的 Lemaitre 应 变 等 价 性 假设 理论 ” , 即 
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o, = 2v0; +[1-D(e)lbkel * a,DCe) 
(19) 
式 中 : v 为 泊 松 比 ; 下 (MPa) 为 初始 弹性 模 量 。 
结合 式 (16) (15) (5) ,可 得 出 本 研究 提出 的 
岩石 统计 损伤 本 构 模 型 , 即 


2vo, + Ege, (e, € £14) 


2vg; 十 Peiexp| 一 pasan] + 


Oy = Eo 


arar e G o odd 
20 


(20) 
2 模型 参数 确定 


三 参数 的 Weibull 概率 分 布 需要 确定 3 个 参数 ， 
其 中 一 个 参数 损伤 国 值 si 可 通过 联 立 Lade-Duncan 
强度 准则 表达 式 和 弹性 胡 克 定律 进行 推导 得 出 , 即 
式 (15) , 剩 下 的 2 个 参数 m 和 eo , 则 是 根据 应 力 -应 
变 曲 线 特征 点 进行 推导 求解 ,其 推导 如 下 。 

当 轴 回应 变 达 到 峰值 应 力 所 对 应 的 轴 向 应 变 
时 ,其 切 向 斜率 为 0, 即 


aa， 


=0 (21) 


ðE: E] F £e 
结合 式 (20) 和 式 (21) 的 微分 结果 化 简 ,得 
mbke. = ma, €. 7 £j ex 

€Q ( ) 


"ef (22) 


E0 
当 s， = e, 时 ,其 应 力 -应 变 应 满足 下 式 。 
gc, = 2v0, * l1 - D(e,)lEe, c a, * Dle.) (23) 
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式 中 , o, MPa) 为 岩石 的 峰值 强度 。 


通过 联 立 式 (22) 和 式 (23) ,得 
Ele, — e, 
em g,-2vo,-c 
(Ee, -ol I (2:03 e) 
€. O | Ee, - 0, | 
(24) 
( E Ze € I yj"! 1/m 
el = | mie, = €, Eld | (25) 


通过 式 (24) 和 (25) 可 知 ,2 个 参数 m 、s 是 关 
于 围 压 os .残余 前 切 强度 or, 峰值 强度 o. 以 及 应 变 
e, 的 函数 表达 式 , 其 中 参数 ou co. 与 o, 存在 明显 
的 非 线 性 关系 ,参数 e. 与 os 存在 明显 的 线性 关 
系 , 即 


O., = aof (26) 
Or 5 Og of" (27 ) 
e. = (ao, + b)/100 (28) 


已 经 有 学 者 对 式 (26) (27) 和 (28 ) 的 适用 性 进 
行 了 相关 的 验证 '。 以 上 3 个 式 子 可 以 很 好 地 描 
述 峰 值 应 力 强度 o., 残余 剪 切 强度 or, 和 峰值 应 力 
强度 所 对 应 的 轴 向 应 变 e. 5 EE os 之 间 的 变化 关 
系 ,并 在 一 定 的 程度 上 简化 了 计算 过 程 。 


3 ”模型 验证 


为 了 验证 本 研究 提出 的 统计 损伤 本 构 模 型 ( 式 
20) 的 适用 性 ,引用 ADACHI! #1 ZHANG 4? x42 
种 不 同 软 岩 的 三 轴 压 缩 试验 结果 与 模型 曲线 进行 对 
比分 析 ,预测 模型 参数 见 表 5 


表 5 预测 模型 参数 


Tab.5 Prediction model parameters 


gos 108 软 岩 2091 
参数 
0.5 MPa 1.0 MPa 1.5 MPa 2.0 MPa 0.5 MPa 1.0 MPa 1.5 MPa 2.0 MPa 
mo 2.138 095 529 2.172577355 1.445730217 0.984961375 1.236412577 1.259212139 1.525 271714 0.634 856 159 
£0 0.010 141095 0.007 908758 0.006940419 0.006 969 095 0.008 889 439 0.009 520286 0.007375 894 0.002219 906 


软 岩 17 的 三 轴 压 缩 试验 偏 应 力 -应 变数 据 和 


着 围 压 施加 的 时 间 增 长 ,岩石 开始 发 生 损 伤 ,其 损伤 


模型 预测 曲线 见 图 4。 可 以 看 出 ,统计 损伤 本 构 模 
型 可 以 较 准 确 地 预测 岩石 损伤 破坏 过 程 及 其 力学 特 
性 ,尤其 是 修正 后 损伤 阐 值 与 实际 试验 数据 也 是 较 
吻合 的 。 在 不 同 围 压 下 的 软 宕 1 损伤 演化 规律 如 图 
5。 从 图 5 中 可 以 得 出 , 宕 石 在 不 同 围 压 下 的 损伤 演 
化 规律 ,在 岩石 处 于 弹性 阶段 时 ,损伤 变量 D =0, 随 


变量 D 也 逐渐 开始 积累 并 最 后 趋 近 于 1, 力 学 性 能 
也 开始 逐渐 衰减 。 在 不 同 的 转 压 下 ,岩石 的 损伤 演 
化 规律 是 有 所 不 同 的 , 随 着 围 压 的 增 大 ,岩石 的 损伤 
闷 值 也 逐渐 增 大 ,岩石 的 损伤 变量 的 变化 速率 逐渐 


软 岩 2 ”的 三 轴 压 缩 试验 的 偏 应 力 -应 变数 据 
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以 及 模型 预测 曲线 见 图 6。 可 以 看 出 ,模型 预测 曲 
线 与 试验 数据 是 基本 吻合 的 , 较 好 地 预测 了 岩石 损 
伤 破坏 的 过 程 。 图 7 则 是 软 岩 2 在 不 同 围 压 下 的 损 
伤 规律 演变 曲线 ,其 损伤 变量 随 轴 办 应 变 逐 渐 增 长 ， 
并 趋 近 于 1。 其 再 一 次 验证 本 研究 提出 的 统计 损伤 
本 构 模 型 的 合理 性 ,并 且 更 进一步 验证 了 其 广泛 的 


me 全 了 试验 值 


8 
模型 预测 

7 

6 
S 
85 围 压 为 20 MPa 
T 099 90,6 4 oo000 
g4 E 291.5 MPa 
iE 


压 为 1.0 MPa 


w 


围 压 为 0.5 MPa 


0 0.02 0.04 0.06 0.08 0.10 
轴 向 应 变 

Kl4 软 岩 1 模型 预测 与 试验 数据 

Fig.4 Soft rock 1 model prediction and test data diagram 

1.0 


器 


0.8 


围 压 为 0.5 MPa 


一 一 围 压 为 1.0 MPa 


一 一 用 压 为 1.5 MPa 


围 压 为 2.0 MPa 


0.2 


0 001 002 003 004 005 0.06 
轴 向 应 变 
图 5 在 不 同 围 压 下 的 软 岩 1 损伤 规律 演化 图 


Fig.5 Damage evolution of soft rock 1 under 


different confining pressures 


9 
m A e 试验 值 
: 模型 预测 


| 732.0 MPa 
Y 


偏 应 力 /MPa 


围 压 为 0.1 MPa 


0 0.02 0.04 0.06 0.08 0.10 
轴 向 应 变 


Ke 软 岩 2 模型 预测 与 试验 数据 


Fig.6 Diagram of soft rock 2 model prediction and test data 
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(0 001 002 003 004 005 0.06 
轴 向 应 变 
7 在 不 同 围 压 下 的 软 岩 2 损伤 规律 演化 图 


Fig.7 Damage evolution of soft rock 2 under 


四 


different confining pressures 


4 结 it 


HA Eb 5 [B [EE IIS 5 T] , V d 5 [d (ELTE TC: 7 4k 
础 ,本 研究 将 Lade-Duncan 强度 准则 与 弹性 胡 克 定 
律 结合 确定 损伤 阔 值 ,对 其 确定 的 损伤 阔 值 的 缺陷 
进行 讨论 分 析 并 进行 了 修正 ,并 根据 修正 后 的 
Lemaitre 应 变 等 价 性 理论 和 Weibull 概率 分 布 理论 建 
立 了 软 岩 的 统计 损伤 本 构 模 型 。 得 出 以 下 结论 。 

1) 本 研究 建立 的 软 岩 统计 损伤 本 构 模 型 ,考虑 
了 损伤 阀 值 对 岩石 力学 特性 的 影响 ,认为 岩石 在 达 
到 损伤 立 值 时 才 会 发 生 损伤 ,其 模型 曲线 较 好 地 吻 
合 了 试验 数据 ,能 较 好 地 描述 软 岩 在 三 轴 压 缩 下 的 
破坏 过 程 和 应 力 - 应 变 关系 ,是 一 个 简洁 且 准 确 的 表 
达 式 。 

2) 通 过 联 立 Lade-Duncan 强度 准则 和 弹性 胡 克 
定律 推导 出 的 损伤 阔 值 求解 表达 式 , 引 用 软 贿 的 三 
轴 压 缩 试验 数 据 对 其 进行 验证 ,其 求解 出 来 的 损伤 
国 值 过 于 保守 。 

3) 针 对 损伤 阔 值 过 于 保守 ,建立 修正 表达 式 对 
其 进行 修正 ,精度 提高 了 21.35% ~37. 85% ,修正 以 
后 的 损伤 赣 值 更 接近 于 实际 。 
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